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SUMMARY 

Three  theories  are  developed  for  computing  the  thermal  stresses 
in  sandwich  cylinders  with  a  continuous  isotropic  core.  These  are  based 
upon  the  following  idealizations:  (1)  dividing  the  cylinder  into  con¬ 
centric  constant  modulus  sub-cylinders  which  are  analyzed  by  thick 
cylinder  theory,  (2)  treating  the  faces  as  membranes  and  the  core  as  a 
thick  cylinder,  and  (3)  using  a  modified  thin-shell  theory  which  accounts 
for  the  variable  modulus.  Results  from  the  theories  are  compared  for  an 
alumina  sandwich  which  is  typical  of  radome  construction.  It  is  found  that 
the  thin-shell  formulation  results  in  a  simple  equation  but  is  only 
accurate  at  large  values  of  the  mean  radius  to  total  thickness  ratio.  The 
membrane-face  equations  agree  veil  with  the  concentric-cylinder  theory 
but  computational  difficulties  are  comparable  to  the  latter  method.  The 
concentric -thick-cylinder  theory  has  therefore  been  selected  to  be  pro¬ 
grammed  for  machine  computation.  It  may  also  be  used  to  analyze  homogeneous 
cylinders  in  which  the  modulus  is  temperature  dependent  and  therefore 
varies  through  the  wall  thickness. 
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INTRODUCTION 

The  aerodynamic  heating  associated  vith  high  speed  flight  produces 
large  temperatures  and  thermal  gradients  in  missile  structures.  Refractor}' 
materials  capable  of  operating  at  high  temperatures  usually  have  low 
strength /weight  ratios  and  are  therefore  structurally  inefficient.  One 
method  of  increasing  their  structural  effectiveness  is  to  employ  them  in 
sandwich  configurations  in  which  the  core  material  is  made  porous  to  de¬ 
crease  its  density.  As  an  example,  laminated  slip  cast  ceramic  radomes 
are  currently  under  development  by  several  companies.  The  electrical  band¬ 
width  and  cost  advantages  that  these  radomes  appear  to  hold  over  homogeneous 
wall  radomes  have  been  discussed  in  Ref.  (1).  The  structural  advantages 
that  such  radomes  have  in  resisting  aerodynamic  and  inertial  loads  are 
also  apparent,  but  the  effectiveness  with  regard  to  thermal  stresses  is 
not  so  obvious  and  requires  quantitative  investigation.  The  work  reported 
herein  was  undertaken  to  develop  the  analytical  methods  to  be  used  in  such 
a  study. 

Weckesser  and  Suess  have  used  the  cylinder  as  an  analytical  model  to 
compare  the  thermal  shock  characteristics  of  various  materials  for  homogeneous 
wall  radomes  (Ref.  2).  They  find  that  the  thermal  stresses  are  critical  in 
the  region  of  transition  from  laminar  to  turbulent  flow.  Since  the  ratio 
of  the  cylinder  radius  to  the  wall  thickness  is  small  at  the  transition 
point  they  use  thick-walled  cylinder  theory  (Ref.  3)  to  compare  stresses. 
Their  results  have  been  found  to  be  in  reasonable  agreement  with  limited 
wind  tunnel  tests  (Ref.  4).  The  cylindrical  configuration  was  also 
selected  for  the  sandwich  investigation  so  that  the  results  could  be  directly 
compared  with  the  previously  determined  data  for  homogeneous  wall  con¬ 
struction. 

The  theories  which  have  been  considered  are  based  upon  the  following 
idealizations  of  the  sandwich  cylinder: 

1.  Dividing  the  cylinder  into  concentric  sub-cylinders  with  constant 
mechanical  properties  in  each  region.  The  regions  are  then 
treated  as  thick-walled  cylinders  and  radial  displacements  and 
stresses  are  matched  at  their  interfaces. 
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2.  The  thin  faces  layers  of  the  sandwich  are  treated  as  membranes 
and  the  core  as  a  thick-walled  cylinder.  Radial  forces  and 
displacements  are  matched  at  the  junction  of  the  faces  and  the 
core. 

3.  The  sandwich  is  analyzed  by  the  Love's  theory  of  thin  shells 
modified  to  include  the  effects  of  the  non-homogeneity  of  the 
wall . 

The  derivations  for  the  above  theories  are  contained  in  Appendices  A 
through  D.  A  comparison  of  the  results  obtained  by  using  the  theories 
is  given  in  the  next  section. 

RESULTS 

In  order  to  obtain  a  comparison  of  the  accuracy  and  relative  com¬ 
putational  difficulties  of  the  three  theories,  thermal  stresses  were 
determined  in  the  sandwich  cylinder  shown  in  Figure  1.  The  geometry  aqd 
mechanical  properties  are  typical  of  an  alumina  sandwich  radome  having  a 
core  with  a  porosity  of  677..  The  temperature  distribution  was  taken  from 
Ref.  (5)  and  is  shown  in  Figure  2.  Mechanical  properties  and  thermal 
expansion  data  for  the  alumina  were  obtained  from  Refs.  (6)  and  (7) 
assuming  that  Poisson's  ratio  and  thermal  expansion  are  not  affected 
by  porosity. 

Circumferential  stresses  (o^)  and  axial  stresses  ( o were  computed 
at  the  inner  radius,  and  radial  stresses  (a^.)  were  found  at  the  juncture 
of  the  inner  face  with  the  core  by  each  of  the  theories.  Results  are 
shown  in  Figure  3  for  a  range  -x  outside  radii  assuming  the  same  thermal 
gradient  and  wall  geometry  in  all  cases. 

DISCUSSION 

The  equations  for  the  concentric  thick-walled  cylinder  are  derived 
in  Appendices  A  and  B.  The  results  in  Figure  3  were  obtained  using  3 
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cylinders,  one  for  each  of  the  face  layers  and  one  for  the  core.  This 
development  of  the  problem  utilizes  the  three-dimensional  theory  of 
elastic  bodies  in  both  the  faces  and  the  core  and  is  therefore  the  most 
exact  of  the  three  theories.  However,  it  is  not  convenient  to  oDtain  a 
closed  form  solution  by  this  method  since  this  would  require  a  general 
solution  of  6  simultaneous  linear  algebraic  equations  tc  evaluate  the 
constants  of  integration.  The  analysis  is  therefore  only  practical  when 
a  digital  computer  is  used  to  solve  the  equations  for  a  particular  design. 

The  concentric  thick-wall  cylinder  solution  is  also  useful  in  the  analysis 
of  homogeneous  cylinders  since  it  is  possible  to  take  into  account  the 
variation  of  modulus  of  elasticity  with  temperature.  This  is  done  by 
dividing  the  homogeneous  wall  into  concentric  cylinders  and  evaluating 
the  modulus  in  each  region  at  the  mean  temperature  for  that  region. 

The  membrane  face  theory  derived  in  Appendix  C  is  restricted  to  those 
cases  in  which  the  face  thicknesses  are  very  small  compared  to  the  radius. 

It  assumes  that  the  radial  stresses  in  the  faces  are  negligible  compared  to 
the  axial  and  circumferential  stresses  so  that  a  two-dimensional  stress- 
strain  law  can  be  used.  It  further  assumes  that  the  axial  and  circumferential 
stresses  are  constant  through  the  face  thicknesses.  Uncertainties  are 
introduced  into  the  analysis  by  lumping  the  faces  into  membranes  at  dis¬ 
crete  radii.  It  is  not  clear  if  the  radial  stresses  and  displacements 
of  the  membrane  and  core  should  be  matched  at  the  outer  surface  of  the 
face  layer,  the  center  of  this  layer,  or  at  its  interface  with  the  core. 

The  results  in  Figure  3j  which  are  for  the  membranes  at  the  interfaces, 
show  good  agreement  with  the  more  exact  concentric  cylinder  theory  over 
the  range  of  radii  investigated.  While  the  membrane -face  theory  results 
in  a  closed  form  solution  the  equations  are  lengthy  and  time  consuming 
to  evaluate. 

The  non-homogeneous  thin  shell  theory  is  derived  in  Appendix  D.  The 
development  assumes  that  normals  to  the  reference  surface  remain  normal, 
straight,  and  unchanged  in  length.  Radial  thermal  strains  are  thereby 
ignored.  Radial  stresses  are  considered  negligible  compared  to  the  axial 
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and  circumferential  stresses  throughout  the  body  so  that  the  stress-strain 
relations  are  two-dimensional.  Equilibrium  is  only  satisfied  on  a  macroscopic 
scale  by  the  stress  resultants  rather  than  on  an  infinitesimal  basis  by  the 
stresses.  The  resulting  closed  form  solution  is  simple  to  evaluate  but  it 
is  seen  from  Figure  3  that  it  is  grossly  in  error  when  the  radius  to 
thickness  ratio  is  small.  For  the  case  which  was  studied  the  error  relative 
to  the  concentric  cylinder  solution  reduces  to  6%  for  a  sandwich  mid-radius 
to  total  thickness  ratio  of  5  and  is  less  than  4 °L  for  a  ratio  of  10.  It  is 
interesting  to  note  from  Appendix  D  that,  except  in  the  edge  regions,  the 
stresses  in  a  thin  shell  of  revolution  do  not  depend  upon  the  shape  of  the 
shell  if  the  temperature  only  varies  through  the  thickness.  Furthermore, 
results  obtained  by  using  an  AVCO  computer  program  (Ref.  8)  show  that  this 
is  also  true  when  the  temperature  variation  in  the  mendional  direction  is 
gradual,  in  which  case  the  stresses  depend  only  upon  the  local  temperature 
distribution  through  the  thickness  and  may  be  computed  from  Eq.  (D.14). 

The  behavior  displayed  by  the  circumferential  and  axial  stresses  in 
Figure  3  at  small  radii  is  interesting.  It  is  seen  that  the  circumferential 
stresses  decrease  but  that  the  axial  stresses  increase.  This  same  situation 
does  not  prevail  in  a  homogeneous  cylinder  where  the  circumferential  and 
axial  stresses  are  equal  regardless  of  the  radius  and  both  increase  with  a 
decrease  in  radius.  The  reduction  of  core  modulus  appears  to  be  effective 
in  decreasing  the  circumferential  stresses  since  it  increases  the  relative 
radial  motion  between  the  inner  and  outer  faces.  However,  significant 
changes  do  not  occur  in  the  relative  axial  motions  except  in  the  edge 
regions  so  that  the  axial  stress  in  the  inner  face  remains  high.  The  lengths 
of  the  edge  region  where  reduced  axial  stresses  occur  will  increase  but 
for  the  length  to  thickness  ratios  used  in  radomes  there  will  still  be 
a  central  portion  where  plane  sections  remain  plane  as  assumed  in  the 
theories.  The  thin  shell  theory  does  not  admit  to  relative  radial  motions 
and  therefore  does  not  show  the  reduction  of  the  circumferential  stresses 
for  small  values  of  the  ratio  of  mean  radius  to  total  thickness. 
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CONCLUDING  REMARKS 

Three  theories  have  been  developed  for  the  analysis  of  sandwich 
cylinders.  As  would  be  expected,  the  modified  thin  shell  theory  results 
in  the  simplest  equation  but  gives  large  errors  at  values  of  the  mean 
radius  to  total  thickness  ratio  less  than  5.  The  membrane  face  theory 
is  in  good  agreement  with  the  concentric  thick  cylinder  analysis  but 
requires  considerable  computational  effort  making  it  impractical  for 
repeated  calculations  unless  it  is  programmed  for  machine  computation, 
in  which  case  the  concentric-cylinder  analysis  is  preferable.  As  a 
result  the  theoretical  analysis  described  in  Appendices  A  and  B,  which  is 
also  useful  in  analyzing  homogeneous  radomes  in  which  the  modulus  is 
temperature  dependent,  has  been  selected  for  machine  programming. 

On  the  basis  of  limited  computations  it  appears  that  the  low  modulus 
core  is  effective  in  reducing  the  circumferential  stress  at  low  values 
of  the  mean  radius  to  total  thickness  ratio  but  that  the  axial  stresses 
are  not  decreased  . 
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APPENDIX  A 

Derivation  of  Equations  for  Concentric  Thick-Walled 
Cylinder  Theory 

The  non-homogeneous  cylinder  is  assumed  to  be  divided  into  concentric 
sub-cylinders  as  shown  in  Figure  4  so  that  the  modulus  of  elasticity  is 
constant  within  each  layer.  To  identify  the  layers  they  are  numbered 
consecutively  from  1  for  the  inner  to  n  for  the  outer  cylinder.  The 
inner  and  outer  radii  of  the  ith  cylinder  are  designated  by  b^  and  b^  _ 
respectively. 

It  is  assumed  that  the  materials  are  isotropic  and  that  Poisson's 
ratio,  v,  is  constant  in  all  layers.  This  assumption  does  not  impose 
severe  restrictions  on  the  theory  for  v  of  most  structural  materials 
(within  the  elastic  limit)  falls  between  0.25  and  0.33.  For  a  given 
material  moderate  changes  in  porosity  do  not  affect  v  and  the  temperature 
dependency  is  slight  so  that  this  assumption  is  fully  justified  in  sand¬ 
wich  cylinders  made  of  a  single  material  but  with  different  porosities  in 
the  various  layers.  Since  E^  is  constant  in  the  ith  cylinder,  we  may 
apply  the  Duhamel  equations  (Ref.  3)  for  thermal  deflections  and  stresses 
in  a  homogeneous  hollow  cylinder  to  the  region.  If  we  assume  the  ends 
of  the  cylinder  to  be  restrained  against  axial  displacement  these  are 
given  by  Eqs.  (c)  through  (f)  from  page  409  of  Ref.  (3).  For  the  ith 
region  we  obtain 
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In  these  equations  and  are  constants  of  integration  which  apply  to 

the  ith  region. 

The  coefficient  of  expansion  has  been  brought  under  the  integral  sign 
in  writing  these  equations  since  it  is  always  permissible  to  consider  aT 
(the  unrestrained  strain  due  to  temperature)  as  a  single  quantity.  As  a 
result  a  need  not  be  linear  over  the  range  of  temperature  changes. 

Since  there  are  n  regions  a  total  of  2n  constants  of  integration  must 
be  found.  These  constants  are  determined  so  that  equilibrium  and  com- 
patability  are  satisfied  on  the  faces  and  interfaces  of  the  layers. 

Equilibrium  on  the  inner  and  outer  faces  requires  that 

ar  0^)  =  -V-l  (A-5) 
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where  p^  and  p  ^  are  the  external  pressures  on  the  inner  and  outer  faces. 

At  each  of  the  interfaces  the  equilibrium  condition 

°r  °W  =  ar  (bi+1);  i=l,2,...,(n^l)  (A-7) 
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and  the  compatability  condition 
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must  be  satisfied.  Substituting  Eqs.  (A-l)  and  (A-2)  into  (A-5)  through 
(A -8)  gives 
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Solution  of  the  2n  simultaneous  equations  generated  by  Eqs.  (A-9)  through 
(A-12)  gives  the  and  coefficients.  Substitution  of  these  into 
Eqs.  (A- 3)  and  (A-4)  yields  aQ  and  a  for  restrained  ends. 
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The  resultant  axial  force  for  restrained  ends  is 
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To  determine  the  axial  stresses  for  unrestrained  ends  we  determine  the 
stresses  due  to  an  axial  force  of  -R  and  superimpose  these  upon  the 
previously  determined  stresses.  The  stresses  due  to  -R  are  (Appendix  B, 
Eq.  (B-l)) . 
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In  the  ith  region  the  radial  and  circumferential  stresses  for  unrestrained 
ends  are  then  given  by  Eqs.  (A-2)  and  (A-3).  The  axial  stress  is 


a 

z , 

l 


aTE.  2vE  ,C . , 

_ i  + _ i-iL-  + 

1-v  (1+v) (l-2v) 


Ei(R/ir) 


ej  (bi+i 


(A-15) 
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The  quantity  defined  by  Eq.  (A-13)  will  usually  have  to  be  evaluated 
by  a  numerical  integration  procedure  such  as  Simpson's  rule.  When  the  wall 


thickness  of  the  ith  cylinder  is  small  it  is  sufficiently  accurate  to  assume 
that  aT  varies  linearly  through  the  thickness  so  that 


aT(r)  = 


jbi+1  aT(bi)  -  b.aT(b.+1)  +[crT(bi+1)  -  aT(bi)]r  (A-16) 


Substituting  this  into  Eq.  (A-13)  we  find  for  a  thin  wall 


Ai=  ( "I(b>)(iaiti-3b2i  ^+i«bl>“ V'2bM-3bibi+i+b3i>]  <A-17) 


This  equation  is  poorly  suited  for  numerical  computations  since  it  results  in 
small  differences  of  large  numbers.  An  equation  of  improved  numerical  accuracy 
may  be  obtained  by  substituting  the  relation  b^+^  =  b^  +  t^  (where  t^  is  the 
thickness  of  the  wall  of  the  ith  cylinder)  into  Eq.  (A-17)  which  upon  simplifica¬ 


tion  becomes 


A.  _ 
1  - 


6(l-v) 


i _ [«T(bi)  (3bi+t.)-krT(bi+1)  (3b.+2t.)] 


(A- 18) 
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APPENDIX  B 


12. 


Derivation  of  Equations  for  Axially  Loaded 
Concentric  Thick-Walled  Cylinders 

In  the  thermal  stress  analysis  for  non -homogeneous  cylinders  given  in 
Appendix  A,  the  stresses  due  to  an  axial  load  were  required.  In  this 
appendix  we  shall  show  by  the  inverse  method  that  Eq.  (A-15)  correctly 
describes  the  stresses. 

The  cylinder,  which  is  the  same  as  described  in  Appendix  A  and  shown 
in  Figure  4,  is  assumed  to  be  subjected  to  normal  forces  distributed  over 
the  ends  which  have  a  resultant  tensile  force  R  directed  parallel  to  the 
axis.  The  cylinder  is  long  so  that,  except  in  the  immediate  vicinity  of  the 
ends,  a  state  of  generalized  plane  strain  exists.  The  stresses  therefore 
only  depend  upon  the  resultant  of  the  end  forces  and  not  upon  their  dis¬ 
tribution.  We  assume  R  to  be  divided  among  the  regions  in  proportion  to 
their  stiffnesses  so  that 


8 

d 


i 

& 


$ 

s 


E.R 

l 


}  (B-l) 


J 

describes  the  stresses  in  the  ith  region.  These  stresses  satisfy  the  equations 
of  equilibrium 


3a 

3t 

ar  '  ae 

+  ri  i=0 
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3r 

3z 

r 

3t 

3a 

T 

rz . 
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+ 

zi 

+  rzi  =  0 

3r 

3z 

r 

which  apply  to  axially  symmetrical 

stress 

distributions  in 

a  solid  of  revolution 

(Eqs.  177,  Ref.  3).  Eqs.  (B-l) 

also  satisfy  the  following 

compatibility 

equations  (Eq.  (g),  page  346,  Ref. 

(3) 

I 

1 


Tha  John*  Hopkin*  Univ«rk»fy 

APPLIED  PHYSICS  LABORATORY 

Silver  Spring,  Maryland 


13 


V2  a  -  1 

ri  r2 


^Gr.  “  +  1+v  sr2 

11 


1  B2® 


=  0 


V2  +  -  (a  -  afi  )  +  -J-  .  - 

9 .  r2  r.  0 .  *  1+v  r  Br 

i  r  i  i 


=  0 


v3  o  2!|  -  o 

z,  i+v  az2 


v2  T  -  -  T  +  — -  —  =  o 

rzi  r2  rzi  1+v  BrBz 


where  ®  is  the  sum  of  the  normal  stresses,  which  from  Eqs.  (B-l)  is  equal 

to  a  ,  and 
z. 

l 


V2 


a!  .  l  a_  ,  af. 

Br2  r  Br  Sz2 


is  the  Laplacian  operator  for  axial  symmetry.  The  boundary  conditions  to 

be  satisfied  on  the  cylindrical  faces  are  given  in  Appendix  A  by  Eqs.  (A-5) 

through  (A-8).  Eqs.  (A-5)  through  (A-7)  are  identically  satisfied  because 

ct  =0.  Since  o  is  the  only  non-zero  stress,  the  strains  are  given  by 
l  zi 

the  uniaxial  stress-strain  laws 


e 


z 


i 


The  strain  displacement  equations  are  (Eqs.  178,  Ref.  3) 


(B-2) 


(B-3) 
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By  simultaneously  solving  Eqs .  (B-l)  through  (B-3)  we  find  the  displacements 
to  be 


Ui  = 


vRr 


n 


TT  £  Mb*  ,-b®) 


j=l 


r  j+i  j' 


wi = 


Rz 


n 


n  L  E.(b*  -b2) 

j=l  J  J+l  J 


(B-4) 


The  first  of  these  equations  satisfies  the  remaining  boundary  condition  (A-8) 
at  each  of  the  interfaces  of  the  cylindrical  surfaces.  The  second  equation 
indicates  that,  as  supposed,  plane  sections  remain  plane.  Finally  we  see 
Eqs.  (B-l)  satisfy  the  equation 


R  =  2tt 


£ 


i=l 


rdr 


thereby  assuring  a  state  of  equilibrium  between  the  applied  force  R  and 
the  internal  stresses 


\ 
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APPENDIX  C 

Derivation  of  Equations  for  Membrane -Face  Theory 

In  sandwich  structures  the  face  layers  are  often  thin  compared  to  the 
core  layer.  This  suggests  a  simplified  theory  in  which  the  faces  are  treated 
as  membranes  while  the  core  is  analyzed  as  a  thick-walled  cylinder.  The 
geometry  and  notation  are  shown  in  Figure  5.  The  following  assumptions  are 
made. 

(1)  The  temperatures  in  the  face  layers  are  uniform  through  their 
thicknesses.  These  temperatures  are  designated  T^  and  T^;  the 
moduli  E^  and  E^  are  determined  at  the  corresponding  temperatures. 
Average  temperatures  are  used  when  there  are  gradients  through 
the  face  layers. 

(2)  The  only  stresses  of  importance  in  the  membrane  faces  are  Qg  and 

a  ,  and  these  stresses  are  constant  through  the  membrane  thicknesses, 
z 

(3)  The  core  region  2  is  a  "thick"  cylinder,  i.e.  the  effects  of  radial 
stresses  are  considered. 

(4)  The  temperature  distribution  in  the  core  is  arbitrary  but  is 
constant  and  is  evaluated  at  the  mean  temperature  of  region  2. 

(5)  The  membrane  stiffnesses  are  lumped  at  a  and  b,  the  inner  and  outer 
radii  of  the  core  respectively. 

As  a  result  of  assumption  (3)  the  deflections  and  stresses  in  the  core 
for  restrained  ends  are  given  by  Eqs.  (c)  through  (f),  page  409  of  Ref.  (3), 
which  are 


1  2 
-  /  aTrdr  +  C,r  +  — 
r  I  1  r 


(C-l) 


0r  “ 1-v 


>  aTrdr  + 


(C-2) 


if-- 
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aTrdr  - 


otTEg  u? 

1-v  1+v 


(C-3) 


Q?TE2  2vE2C1 

°z  ~  1-v  +  (1+v) (l-2v) 


(C-4) 


From  Eqs .  (C-l)  and  (C-2),  the  radial  displacements  and  stresses  at  a  and 
b  are  given  by 


u  (a)  =  a  +  _2 


(C-5) 


u  (b)  =  (1+v)  ^  +  ci  b  +  ~ 


(C-6) 


a.  (a)  =  hi  A. 
1+v  l-2v 


(C-7) 


CT  (b)  V  A 

O  '  i  .  .  . 


1+v  \  l-2v  b2 


(C-8) 


where 


A  =  T^J  /  aTrdi: 


(C-9) 


Because  of  assumption  (2)  the  strains  in  the  inner  face  are  given  by 
the  two-dimensional  stress-strain  equations 


‘ei  'i  (°0i  •  V°2J+  (o,T)i 


li’il  [°zi  -  v°9i)+ 


(C-10) 
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where  v  is  assumed  to  be  the  same  as  for  the  core.  With  restrained  ends 
=  0  and  it  follows  from  the  second  of  Eqs.  (C-10)  that 

°Z1  =  V\  "  (CVT)1  E1  (C“U) 
Substituting  this  into  the  first  of  Eq.  (C-10)  gives 


\  =  ^E^  X  +  <1+V^aT>l  (G“12) 

The  circumferential  strain  in  the  inner  membrane  face  is  related  to  its 
radial  displacement  u^  by 

e91  =  ^1  (C-13) 

a 

The  radial  stress  in  the  core  at  r  =  a  imposes  a  bursting  pressure  on  the 

inner  membrane  so  that  ctq  is  given  by  the  hoop  stress  equation 

91 

a.  .  £E  .  k(a)  *  Pl1 a 

91  t  tx 

where  p^  is  the  inside  pressure. 

From  Eq.  (C-7)  this  may  be  written 

G1  °2|  .  Pla  (C-14) 

l-2v  a2/ 

£2.1  (C-15) 

b 3 1~  H 

where  p^  is  the  outside  pressure. 


In  a  similar  fashion,  for  the  outer 

face  we 

find 

F,„A 

E0b 

C, 

gn  2 

3'bt3  ' 
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t3(l+v) 

l-2v 

E„a 

CTe  =  2  . 

1  t, (l+v) 
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By  combining  Eqs .  (C-12)  through  (C-14)  we  obtiiin 


18. 


ui  a 


v_ 

E1t1 


(1-v) 


W1 

l-2v 


A  +  Et 


Pla 


(1-v3)  +  (1+v)  (aT) 


11 


•1 


(C-16) 


and  by  the  same  procedure  the  radial  displacement  of  the  outer  face  is 


u3  =  b 


(1-V3) 


V 

E3t3b 


(1-v) 


v_ 

E3t3 


“1 

l-2v 


2 

l-2v 


Pob 

(l-vs)+(l+v)(aT)3 


(C— 17) 


Noting  that  u^  =  u(b)  for  compatible  deformations  of  the  faces  and  core  we 
find  the  following  simultaneous  equations  for  and  by  equating  Eqs. 
(C-5)  with  (C-16)  and  (C-6)  with  (C-17) 


DC1  +  EC2  =  J 
FC1  +  GC2  -  K 


where 
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Solving  the  simultaneous  equations  leads  to  the  results 


19. 


=  (JG  -  EK)  / (DG  -  EF) 
C2  =  (DK  -  FJ)/(DG  -  EF) 


(c- 


18) 


The  a  stresses  in  the  core  for  restrained  ends  are  obtained  by  sub- 
z  J 

stituting  into  Eq.  (C-4).  The  resultant  restraining  force  in  the  core 
is  found  by  integrating  Eq.  (C-4)  over  the  core  area  so  that 


rdr 

z 

vE2C1(b2-a2)' 
(1+v)  (l-2v) 


For  restrained  ends,  the  axial  stress  in  the  inner  face  is  found  by  sub¬ 
stituting  Eq.  (C-14)  into  Eq.  (C-ll)  giving 


vaE„ 

7  _  2 

Z1  (l+v)t, 


~1 

l-2v 


(aT)1  Ex 


(C-19) 


which  produces  a  resultant  force 


=  2nat 
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vE2a 


(l+v)t..l  l-2v 


-  -4  -  (^t) 
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In  a  similar  fashion  we  find 


Z3  = 


vE, 
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1+v  \l-2v  b 
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The  total  restraining  force  is  then 

R  =  .Rx  +  +  R3 


which  gives 


J  =  "2  [(1-v)E2A  +  (aT)1E1at1+(aT)3E3bt3]  (C-21) 

The  stresses  for  unrestrained  ends  are  determined  by  superimposing  the 
stresses  for  an  axial  load  of  -R.  By  following  the  method  of  Appendix  B 
we  find  these  stresses  to  be 


Ei(R/rr) 

°z  =  2E]at]+E2 (b* -a* )+2E3bt3 


>  (C-22) 


o 

r 


ae 


=  o 


) 


The  total  stresses  for  unrestrained  ends  are  obtained  by  substituting 
che  stresses  from  Eqs.  (022)  upon  the  previously  determined  stresses  for 
restrained  ends.  The  results  are:  for  the  inner  face 


a 

zi 


vE9a 

U+vTt, 


l-2v 


(l+v)t1  \  l-2v 


-  (»T)1E1 


E1(R/tt) 

2E1at1+E2  (ba -a*  )+2E3bt3 


(C-23) 
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Derivation  of  Equations  for  Nonhomogeneous 
Shell  Theory 


I 

m 


In  this  appendix  the  equations  for  the  stresses  and  displacements  in 
a  nonhomogeneous  thin  shell  of  revolution  subjected  to  a  temperature  dis¬ 
tribution  which  only  varies  through  the  thickness  are  derived.  The  non¬ 
homogeneity  is  assumed  to  occur  in  the  thickness  direction  only. 

The  thermal  stresses  give  rise  to  the  stress  resultants  shown  in 
Figure  6.  These  are  defined  as  follows 


N  =  la 

cp  J  <p 


M  =  fa 
cp  J  cp 


dz 


zdcp 


W 


“e  '  J  °t 


dz 


Me  = 


zdz 


T .  dz 
i|fz 


)  (D.l) 


where  cp  and  9  are  coordinates  shown  in  Figure  7  and  z  is  a  coordinate 
measured  from,  and  normal  to  the  reference  surface  (defined  later)  of  the 
shell.  The  integrals  extend  over  the  shell  thickness  h. 

Force  and  moment  equilibrium  conditions  for  the  differential  element 
of  Figure  6  lead  to  the  following  equations  (pp.  534,  Ref.  9) 


(N  r  )  -  N  r  cos  cp  -  r  Q,.  +  r  r. 
dtp  tp  o  01  o  cp  ol 


Y  =  0 


d(Q  r  ) 
vxCPo 


N  r  +  NQ  r.  sin  cp  +  — t*  ■  +  Zr.r  =  0 

cp  o  0  1  T  dtp  lo 


—  (M  r  )  ■  M0  r.  cos  cp  -  Q  r.r  =  0 

dtp  '  tp  o;  0  1  Y  xcp  1  o 
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where  the  radii  rQ  and  r^  are  shown  in  Figure  7.  No  external  loading  is 
applied  so  that  the  surface  forces  Y  and  Z  in  these  equations  are  zero. 

The  material  is  considered  to  be  linearly  elastic.  In  addition  the 

following  assumptions  which  are  usual  in  shell  theory  are  made:  (a)  the 

material  is  isotropic  in  the  surface  of  the  shell,  (b)  the  modulus  of 

elasticity  in  the  z  direction  and  the  shear  modulus  which  relates  t  to 

zqj 

are  infinite,  and  (c)  the  coefficient  of  expansion  in  the  z  direction, 
and  the  Poisson's  ratio  which  relates  the  strain  in  the  z  direction  to  the 
strains  in  the  surface  directions  of  the  shell  are  zero.  This  is  equivalent 
to  the  Kirchoff  assumptions  that  normals  to  the  reference  surface  remain 
normal  and  unchanged  in  length  during  deformation  of  the  sh6ll.  As  a 
result  the  plane-stress  Hooke's  law  equations 


%  *  i  <V  w8)  +  “T 


Ee  *  i  <ae  -  %>  +  “T 


apply.  Solving  these  for  the  stresses  gives 


(D.3) 


E  N  o>ET 

o  =  - — 2  (e  +  vea)  -  - — 
cp  1-v  cp  9  1-v 


_  E  ,  .  N  aET 

0  1-v  N  6  9  1-v 


(D.4) 


Using  the  previously  mentioned  shell  assumptions  that  normals  to  the 
reference  surface  remain  normal  and  unchanged  in  length  during  deformation 
of  the  shell,  the  strains  may  be  expressed  as  (p.  431,  Ref.  9) 


e  =  e  (0)  -  x  z 
cp  cpW  Acp 


e6  =  «6<0)  -  Xez 


(D.5) 


where  e  (0)  and  eo(0)  are  the  strains  in  the  reference  surface  (z=0)  and 

CD  u 

X .  and  Xq  are  the  changes  of  curvature  of  the  reference  surface. 

Cp  W 
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Placing  these  into  Eq.  (D.4)  gives 


°«p  ■  tV  [  y°> +  «8<®>  - 2  % +  vl  - 

ce  =  [  V°>  +  'V°)  '  2  &e  +  %>] 


By  substituting  the  equation  for  into 
from  Eqs.  (D.l)  we  find 


the  relationship  for  N  and  M 

cp  cp 


%  +  vxe) 


%  +  v><0) 


XD.6) 

Zy 


h 


where  it  is  assumed  that  E  may  vary  through  the  shell  thickness.  These 
equations  are  simplified  if  z  is  measured  from  a  surface  such  that 

/zEdz  =  0.  The  distance  to  this  surface  from  an  arbitrary  initial  surface 


is  given  by 


(D.7) 


where  z  is  the  distance  measured  from  the  arbitrary  initial  surface  and  E, 
o  1 

is  an  arbitrarily  chosen  reference  modulus.  When  E  is  piecewise  constant 
within  layers  of  the  shell  the  last  equation  may  be  written 
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i 


(D.8) 
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In  this  E^  is  the  modulus  in  the  ith  layer  which  has  a  thickness  tu.  The 

distance  z  is  measured  to  the  middle  surface  of  the  ith  layer  and  the 
°i 

sums  extend  over  all  layers.  We  note  that  z'  is  the  distance  to  the 
modulus -weighted  centroidal  surface.  When  the  shell  is  homogeneous  z' 
reduces  to  z,  the  distance  to  the  middle  surface  of  the  shell. 

Measuring  z  from  the  reference  surface  Eqs.  (D.6)  become 

vk’  [y0)+'V0)]  -Pv  (D-9a) 

Mr 

Mv  -  -D’  %  +  *e>  -  TTv  <D-9b> 


and  in  a  similar  fashion 


Ne  =  K'U(0)  +  vV0)]'^v 


Mi 

mq  =  -D'  (Xq  +  VX„)  -  rr. 


e  V  i-v 


(D.9c) 

(D.9d) 


where  the  following  definitions  apply 


K'  El  fE  H 
K  =T^F  Jij  dz> 

h 

Nt  =  J  a  ETdz, 


>  I, 


dz 


mt  =  J 


?  (D.10) 


OfETzdz 


The  quantities  K'  and  D'  are  the  extensional  and  bending  rigidities  of  the 
nonhomogeneous  shell.  When  E  is  piecewise  constant  in  layers  the  first 
two  of  Eqs.  (D.10)  may  be  written 
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k'  -  -X  T  -i 

jSl 

E,  n  E, 


i=l  1 


hi3 
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When  the  shell  is  homogeneous,  K'  and  D’  reduce  to  F  a  ad  D,  the  usual 
shell  stiffnesses  defined  by  (Ref.  9) 


K' 


D* 


=  D  = 


Eh3 

12(1-0 


The  displacement  of  an  arbitrary  point  in  the  reference  surface  is 
resolved  into  two  components:  v  which  is  tangent  to  the  meridian  through 
the  point  and  w  in  the  direction  of  the  inward  normal  to  the  reference 
surface  (Figure  8).  The  strains  and  changes  of  curvature  of  the  reference 
surface  may  be  expressed  in  terms  of  the  displacements  as  follows  (pp. 

534  -  535,  Ref.  (9) 


e 

<P 


(0) 


1  dv  w 

r1  dcp  rx* 


eQ(0)  =  l  cotcp  -  2 
2  2 


-  A  d_  Z  4.  dw 

Xcp  rx  dcp  vl  rjdrp 


N 

)  (D.12) 


where  is  the  principal  radius  of  curvature  shown  in  Figure  7.  Substituting 
Eqs.  (D.12)  into  (D.9)  we  find 


By  substituting  these  into  Eqs.  (D.2)  the  equilibrium  equations  may  be 
expressed  in  terms  of  v,  w,  and  Q  .  A  solution  to  these  equations  is  then 
the  solution  to  the  problem. 
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The  solution  for  the  stresses  in  a  homogeneous  clamped  flat  plate 
of  uniform  thickness  subjected  to  a  temperature  distribution  which  only 
varies  through  the  thickness  is  the  same  as  for  a  cylinder  or  sphere  with 
the  Same  thickness  and  temperature  distribution  (Ref.  9).  One  would 
therefore  suspect  that  the  stress  would  be  the  same  for  any  shell  of 
revolution  and  that  a  similar  situation  may  also  be  true  for  nonhomogeneous 
plates  and  shells. 

The  thermal  stresses  in  a  nonhomogeneous  plate  with  clamped  edges  and 
constant  K' ,  D1,  N^,,  and  M^,  is  (Ref.  10) 


CT 

x 


CT 


y 


E 

l-v 


T  =  T  =  0 
yz  zx 


One  would  be  led  to  believe  that 


CT 


<P 


=  JL  1 

0  i-v  l-v* 


(D.14) 


may  be  the  solution  for  a  nonhomogeneous  shell  of  revolution  in  which  K' , 
D',  N^,  and  are  constant.  It  will  be  shown  that  this  is  true. 


The  stress  resultants  associated  with  Eq.  (D.14)  are  found  by  sub¬ 
stituting  them  into  Eqs.  (D.l),  to  give 


N  =  NQ  =  0 
cp  0 


M 

<P 


l-v 


(D.15) 


=  0 


These  identically  satisfy  the  first  two  of  Eqs.  (D.2)  and  the  third  reduces 
to 


**T  ^ro  r.  coscpj  =  0 

l-v  dcp  “  1  J 
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From  Figure  9  we  see  that  if  there  is  no  discontinuity  in  the  tangent  to 
the  meridian  curve  dro  =  r^cosep  dtp,  or 

dr 

_ o  _  r.  cos  tp  =  0 

dep 

so  that  the  last  equation  is  also  identically  satisfied  and  Eq.  (D.14)  is 
the  solution. 

The  deflections  are  found  by  integrating  the  strains.  By  substituting 
Eq.  (D.14)  into  (D.3)  we  find 

1  Nt 

%  =  ee  =  TZ?  •  F 

which,  from  the  first  2  of  Eqs.  (D.12)  gives 

N„ 


1  dv  w  _  1  “T 

r1  dtp  "  r1  l-v^  *  K' 


v  w  1 

—  cot  CP - -  q — S 

r2  r2 


h 

K' 


(D.16) 


By  eliminating  w  from  these  equations  we  obtain  the  following  differential 
equation  for  v 


dv  1  JT  /  \ 

-  V  cot  <p  -  r  <rl  -  r2> 


dtp 

Integrating  (pp.  446  -  447,  Ref.  9)  we  find 
v  -  sin  cp 


(D.17) 


i  f(i 

T^?  *  K'  J 


r2)  dtp 


sincp 


+  C 


] 


where  C  is  a  constant  of  integration  to  be  determined  from  the  known  value 
of  v  at  the  support.  By  substituting  Eq.  (D.17)  into  the  second  of  Eqs. 
(D.16)  the  following  equation  for  w  is  obtained 


w  =  i^6,F{cos»[/<1si4>9  +cJ‘r2] 


(D.18) 
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Fig.  2  TEMPERATURE  DISTRIBUTION 
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Fig.  6  GEOMETRY  AND  STRESS -RESULTANTS  IN  A  THIN  SHELL 
OF  REVOLUTION. 
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